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1. Introduction 

For a long time, theoretical studies of atomic transition probabilities have mostly 
concentrated on electric dipole (El) allowed transitions, usually responsible for the 
strongest lines in atomic spectra. It has been quickly realized that forbidden transitions 
such as electric quadrupole (E2) and magnetic dipole (Ml) transitions would gain 
in intensity in astrophysical and low-density laboratory plasmas where collisional 
deexcitation of metastable states is low enough to buildup the population of metastable 
states [21 [3]. These forbidden lines are indeed observed in a wide range of astronomical 
objects such as planetary nebulae, WolfRayet stars, novae and the Sun, and they are 
important in providing information on the electron temperature and density of stellar 
objects. In some situations, these lines can be used for calibrating spectrometers over 
wide wavelength ranges [Ij. It is now recognized that more than half of the atoms in the 
Universe recombined via forbidden channels, so that their accurate treatment is crucial 
in order to follow the cosmological recombination process with the level of precision 
required by future microwave anisotropy experiments [5]. 

Amongst the forbidden processes, magnetic dipole transitions often played a key 
role. A review on Ml transitions in Hydrogen and H-like ions has been written 
by Sucher |6]. Helium-like systems occupy a special position in atomic physics as 
the simplest multielectron systems for testing the theoretical calculation of transition 
amplitudes [7j. As illustrated by Johnson et al. in their review [8], a great deal of 
attention has also been given to the magnetic dipole transitions 2^5*1 — 1^5*0 along the 
helium isoelectronic sequence in the period 1970-1995. Historically, the observation of 
the spectral line associated to the transition ls2s — > Is^ ^5* in helium-like ions was 
explained by Breit and Teller [9] as a two-photon electric dipole (2E1) process, the Ml 
transition amplitude being strictly zero when using the "classical" transition operator. 
Later, Gabriel and Jordan [TT] identified some coronal solar lines as arising from the 
transition ls2s — > Is^ in hehum-like C V, O VII, Ne IX, Mg XI and Si XIII, and 
proposed that ls2s would decay to the ground state mostly through the Ml channel. 
Griem [12] demonstrated theoretically on the basis of relativistic calculations that this 
conjecture was founded. 

But the interest of magnetic dipole transitions is not limited to one- and two- 
electron systems. After Bowen's interpretation of the strongest nebular transitions |13j . 
Condon [14J found that the nebulium lines Ni and N2 of O HI were essentially due 
to magnetic dipole radiation. The theoretical tools for estimating Ml transition 
probabilities were then developped by Pasternack [TSj and Shortley [16] in LS coupling 
and in the intermediate coupling. The solar coronal lines remained a longstanding 
spectrocopy mystery until Edlen [T^ found the wave numbers coincidences with Ml 
transitions within the ground configurations 3s^3p, 3s^3j9^, 35^3]?"^ and 3s^3p^ in Fe and 
Ni [18]. A bit later. Ml transitions were identified by Edlen [19] in the configurations 

and and Qp^ of Xenon and Radon. Forbidden lines have now long been used 
in the analysis of solar, stellar or nebular astrophysical plasmas. In the infrared, they 
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are amongst the strongest stellar lines, and many of those of the light atoms have been 
measured. Spectra of the solar corona and solar flares are rich in Ml lines in the visible 
and ultraviolet for the elements nitrogen through nickel. From the eighties, these lines 
have achieved new importance for diagnostics of low-density, magnetically-confined, 
high-temperature plasmas generated in tokamaks, as illustrated by the compilation 
of Kaufman and Sugar [20] of the forbidden lines within ns'^np'^ configurations. The 
majority of the reported observed lines belonging to elements with Z < 20 were observed 
in astrophysical spectra, while the lines belonging to elements with Z > 28 were 
observed in spectra of low density, magnetically confined, high temperature laboratory 
plasmas. This compilation has been extended recently by Feldman and Doschek [21] 
including lines in the 500 — 1600 A range emitted by solar coronal plasmas and recorded 
by the Solar Ultraviolet Measurements of Emitted Radiation (SUMER) spectrometer. 
Magnetic dipole transitions within the lowest 3d^ term in the Ti I isoelectronic 
sequence, violating the general wavelength scaling behavior (ie moving rapidly to shorter 
wavelengths as Z increases) and therefore suitable for plasma diagnostics in high-energy 
Tokamak fusion devices, have been found by Serpa [22], [231 l2l|- 

The theory of magnetic dipole transitions in crystals has been often overlooked in 
the literature but their importance in the description of the spectroscopic properties of 
lanthanides in crystals has been stressed recently [25] . 

In the next two subsections of the introduction, we illustrate the increasing interest 
in M1/E2 transitions of both experimentalist and theoretician communities by providing 
a list, far to be exhaustive, of relevant publications that appeared during the last few 
years. A commented bibliography for forbidden transitions, including Ml transitions, in 
atoms and ions can be found in Biemont and Zeippen [26] covering the period 1989-1995, 
a helpful complement to their previous bibliography on the subject |27j . 

1.1. Experimental work on Ml transitions 

After the theoretical review work by Johnson et al [8], experimental data became 
available for various low-Z elements [28], [29], [30], as illustrated by Trabert et al. [S] 
who explored the heavy-ion storage ring technique for investigating the metastable level 
ls2s ^Si of Li+ and Be^"^. Radiative lifetimes in the microsecond regime were measured 
using the Electron Beam Ion Trap (EBIT) facility, allowing the precise measurement of 
the lifetime of the ls2s ^Si metastable level in heliumlike 0"^+ [28], N^+, 0^+ [29] and 
S^^"*" [32]. This kind of experimental work is crucial for testing the atomic physics theories 
and is not limited to two-electron systems. These experiments provide important checks 
on the theoretical calculations of the rate of the transitions and are not limited to 
heliumlike systems. The precision of measurements of M1/E2 transitions between the 
fine-structure levels in highly charged ions has been continously increased. The lifetimes 
of the P levels producing the coronal transitions of Fe XIV (Al-like) and Fe X (Cl-like) 
have been measured in an electrostatic Kingdon trap [33], producing the metastable 
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ions with an electron cyclotron resonance ion source. The same technique [31] has 
been used to obtain the magnetic dipole transition rates from measured lifetimes of 
levels of Be-like and B-like ions [35], and to measure the lifetimes of metastable levels 
in the ground term of Fe ions within Ss^Sp'^, {k = 1 — 5) [36]. Excited-configuration 
metastable level lifetimes have also been obtained in Cl-like Mn IX and Fe X [37|. A 
similar experimental setup has been adopted for the measurement of the metastable 
lifetime of 2s^2p^ ^5*0 in O^^ [38j. The electron beam ion trap technique [39] and/or 
the heavy-ion storage ring measurements [lOl HI] have been explored to investigate the 
El-forbidden transition probabilities between the fine-structure levels of the ground 
configuration or to measure the El-forbidden excited levels decay rates in Be-like 
[121 SSI E], B-hke [13 sal m SHI mi Sa SI], C-Uke [Sni ETJ, N-Uke [SU], F-like 
[431 SU, Al-like [521 [531 El [55], Si-like [45l IMllSTl [58l [56l [55] , P-hke [58l|59l [60l|55], 
S-hke [nHIMl ESI [601 [55], Cl-like [581 [631 |53ll59], Ti-like [231 El] ions. The accuracy 
achieved in the EBIT measurement of the ls^2s^2p ^-Pg/a level lifetime in boronlike 
Argon [IH] is such that it becomes sensitive to the relativistic electron correlation like 
the frequency-dependent Breit interaction, as well as the electron anomalous magnetic 
moment. From the lifetimes of 35^3^^^ ground configuration levels of A1-, Si-, P- and S- 
like ions of Fe, Co and Ni measured at a heavy-ion storage ring, Trabert et al [55] 
observed that the decay curves show strong evidence of cascade repopulation from 
specific 3d levels, implying that the clean monoexponential (or few-exponential) decay 
signals observed in a number of earlier measurements at heavy-ion storage rings are not 
generally encountered. Atomic level lifetimes have been measured for magnetic dipole 
transitions in Si-like, P-like, Ar-like, K-like, and Ca-like ions of Kr, using an electron- 
beam ion trap in magnetic trapping mode [65] . The storage ring technique has been used 
for measuring radiative lifetimes of metastable levels in more complex atomic systems 
such as Eu II [66], Fe II [671 [68], La II [69], Ti II [70], Ar II [71j, Sc II [72j , for which 
the selectivity in populating the closely spaced metastable levels is achieved using the 
laser probing technique [73] . 

1.2. Theoretical calculations 

Einstein rate coefficients have been calculated for magnetic dipole, electric quadrupole 
and magnetic quadrupole transitions in Be- like ions [71] , using configuration interaction 
wavefunctions in the Breit-Pauli approximation, completing previous studies covering 
selected ions [751 US]- A simple expression for calculating the line strengths of the 
magnetic dipole transitions has been derived by Kingston and Hibbert [77] and tested 
in Be-like ions. More recently, Breit-Pauli energy levels, lifetimes, and transition 
probabilities obtained from variational MCHF wave functions, have been compiled 
for the beryllium-like to neon-like sequences [78]. Breit-Pauli transition probabilities 
were computed by Froese Fischer and Rubin [79] for the intraconfiguration E2 and Ml 
transitions in Fe"*"^, an ion that is expected to be a significant fraction of gaseous iron 
in many nebulae. Relativistic energy levels, lifetimes and transition probabilities have 



Irreducible tensor-form of the relativistic corrections to the Ml transition operator 5 

also been obtained for the sodium-like to argon-like sequences [Sn|, also using the full 
relativistic multiconfiguration Dirac-Hartree-Fock (MCDHF) method to ascertain the 
reliability of the Breit-Pauli approximation at high-degrees of ionization. Radiative 
rates for electric quadrupole and magnetic dipole transitions among the lowest 32 levels 
belonging to the 3d^ configuration of doubly charged iron (Fe III) were calculated by Deb 
and Hibbert |8T] using J-dependent configuration interaction wavefunctions calculated 
within CIVS. Ion storage ring lifetimes measurements of metastable levels in Sc II were 
assisted by Breit-Pauli CI calculations [72J. 

The many-body perturbation theory (MBPT) was used to account for the electron 
correlation, including relativistic corrections in the Breit-Pauli approximation for 
evaluating E2 and Ml transition rates along the N I isoelectronic sequence |82]. As far as 
Ml transition amplitudes are concerned, for all these Breit-Pauli calculations whatever 
the method used for obtaining the zero-order non-relativistic wave functions (variational 
MCHF, MBPT, CI, ...), the relativistic corrections to the Ml transition operator 
that should be taken into account according to Drake [1], have been systematically 
omitted. There is a large number of theoretical works on Ml transition probabilities 
using SUPERSTRUCTURE [83] and its extension AUTOSTRUCTURE [8ll|85], 
also adopting the Breit-Pauli approximation but including the relativistic corrections in 
the transition operator, after the seminal study of Zeippen [86] on the subject. Some 
contributions are presented in section 2, dedicated to that issue. 

In the full relativistic scheme, the relativistic many-body perturbation theory 
(MBPT), including the Breit interaction, was applied to estimate transition rates for 
all Ml transitions within 2/2/' and for some 2/3/' — 2/3/' transitions in Be-like ions 
with nuclear charges ranging from Z = 4 to Z = 100 [87\. Relativistic MBPT 
calculations have also been performed for multipole (El, Ml, E2, M2, E3, and M3) 
transition wavelengths and rates between 3/~^4/' excited and ground states in nickel- 
like ions [88]. This theory agrees with MCDHF calculations in lowest-order, includes 
all second-order correlation corrections, and includes corrections from negative energy 
states. The 4s^4p ^Pi/2-3/2 transition rates were calculated using a similar method for 
Ga-like ions with nuclear charge Z ranging from 31 to 100 [89]. All these calculations 
provide theoretical benchmarks for comparison with experiment and theory. Transition 
rates and line strengths are calculated for electric-multipole (E2 and E3) and magnetic- 
multipole (M1,M2 and M3) transitions between 35^3/3^^°, 3s^3p^3d^4l, 3s^3p^3d^°Al, 
and 3s3p^3d^^4l states with 4/ = 4s, 4p, Ad and 4/ in Ni-like ions with nuclear charges 
ranging from Z = 34 — 100 [90] . More recently, the same approach was adopted to study 
the Ml contributions to the lifetimes of the 5C/3/2 and 5^5/2 metastable levels of Ba"*" 
[9T] . Motivated by the need of accurate theoretical and laboratory data on high-energy 
radiative emission of highly ionized tungsten ions in EUV and Xray wavelength regions, 
Vilkas et al |92j used the relativistic multireference M0ller-Plesset perturbation theory 
for calculating energy levels and transition probabilities in Ne-like xenon, tungsten, and 
uranium ions. 

The variational multi-configuration Dirac-Fock method has been used to study the 
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2s^2p ^-P3/2-i/2 Ml transitions of B-like and all the Ml transitions of Be-like argon |93j . 
On the basis of a multiconfiguration Dirac-Fock method, a systematic study has been 
carried out for the decay process of the ls2s^ ^>S'i/2 state of Li-like ions to discuss 
the balance between the Auger decay channel, the two-electron one-photon (TEOP) 
and Ml radiative decay modes [9l]. Radiative rates are calculated for E1,M1,E2,M2 
transitions in Cl-like Fe X [95] and in Boron- like to F-like Kr ions [96j, using the GRASP 
code. A multiconfiguration Dirac-Fock method has been applied [97J for computing the 
wavelengths and transition probabilities for lines in the X-ray spectra of the gallium-like 
ions from Yb XL to U LXII, including the forbidden transitions (Ml and E2) within the 
ground configuration As^Ap. The E2 and Ml transition probabilities within the ground 
configuration 3(i^ of Fe^"*" have been calculated by Froese Fischer et al |98] , adopting the 
Multiconfiguration Dirac-Hartree-Fock approach. 

The magnetic-dipole transition probabilities between the fine-structure levels for 
B-like and Be-like ions were calculated by Tupitsyn et al. [99], using the configuration 
interaction method in a Dirac-Fock-Sturm basis, getting a fair agreement with 
observation, except for B-like Ar^'^+ for which the discrepancies cannot be explained 
by the recoil corrections [lOOj . 

2. Relativistic corrections to the Ml operator 

From a theoretical point of view. Ml amplitudes played a crucial role in the development 
of theoretical atomic physics. For example, the role of the negative continuum has 
been illustrated in the full relativistic calculation of Ml transitions in two-electron ions 
using the Multiconfiguration Dirac-Fock Method jlOlj . Large contributions of negative- 
energy states to forbidden magnetic-dipole transition amplitudes have been also found in 
relativistic MBPT calculations in beryllium-like ions [87] and alkali-metal atoms |102j . 
The estimation of Ml transition amplitudes in the Breit-Pauli approximation also raises 
interesting questions. Starting from the Dirac-Breit Hamiltonian and the semi-classical 
radiation theory, Drake [I] derived the expression of the relativistic magnetic dipole 
transition operator to be used in the Breit-Pauli scheme. He estimated the contribution 
of this decay mode to the radiative lifetime of the metastable ls2s ^5* level along the 
helium isoelectronic series in the range Z = 2 — 26, in good agreement with observation. 
The relativistic corrections to the magnetic dipole moment operator in the Pauli 
approximation were derived independently by Drake [1] using semiclassical radiation 
theory and by Feinberg and Sucher |103] using conventional quantum electrodynamics. 
In their work, the last authors criticized Drake's approach for including both the 
Breit operator and transverse photon in the same Hamiltonian. In response to these 
criticisms, Drake |104j derived theorems showing that the Breit interaction remains valid 
in the presence of radiation emission and that the 0(«^Z^) corrections to radiative- 
transition probabilities obtained by the semiclassical method always agree with the 
quantum-electrodynamic results. A few years later, Lin [105] obtained a Foldy- 
Woulthuysen transformation for use in the field theory of quantum electrodynamics 
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and clarified the remaining ambiguity concerning the treatment of the term from 
the interaction Hamiltonian. A completely different approach based on relativistic 
many-body perturbation theory was adopted by Johnson et al. [8] and Derevianko et 
al. [106j . Some disagreement with nonrelativistic results based on the Breit Hamiltonian 
incited Lach and Pachucki [107j to derive rigorously nonrelativistically forbidden single- 
photon transition rates between low-lying states of the helium atom within quantum 
electrodynamics. More recently, Pachucki [108j extended this study to El, Ml, E2 and 
M2 forbidden transitions in light atoms, also investigating the role of the anomalous 
magnetic moment of the electron in some forbidden transitions [||. 

The relativistic corrections to the Ml transition operator derived by Drake [1] were 
shown to be relevant for more complex systems of astrophysical interest. The ratio of 
line intensities of forbidden transitions ^-^5/2 ~ ^-^3/2 ~ ^"^3/2 within the half- 

filled atomic configuration provides an interesting tool for plasma diagnostic (see for 
example |110l 1111] ). In the high electron density limit, this ratio is entirely determined 
by the radiative transition probabilities and given by 

To explain the disagreement between the theoretical and observed values of this 
ratio in highly dense plasmas for atomic systems such as O II |112j and N I 
[113] . Zeippen |114j suggested the importance of the relativistic corrections to the 
magnetic dipole transition operator, as originally derived by Drake [1] for two-electron 
systems. Introducing these corrections in the code SUPERSTRUCTURE, Eissner and 
Zeippen [83j and Zeippen [86] calculated the ratio [O II] /(3729)//(3726). This value 
is considered as the best ratio value to be used for estimating electron densities from 
observational spectroscopy [115] . Further calculations have been performed by Zeippen 
and collaborators for O II j ll6] and N-like ions |117l 1118] , getting improved radiative 
transition probabilities for forbidden lines in closer agreement with MCHF/Breit-Pauli 
values [119] but systematically illustrating the effect of the relativistic corrections to 
the Ml operator. The low- and high-density limits of the forbidden fine structure 
line intensity ratio [O II] /(3729)//(3726) have been reinvestigated by Pradhan et al 
|120] using the Breit-Pauli R-matrix (BPRM) method, in considerable disagreement 
with the elaborate calculations of McLaughlin and Bell |121] that Keenan et al (1999) 
|122] used to get a problematic low-density limit of 1.5. Montenegro et al |123j showed 
that this outstanding discrepancy between the observed line intensity ratios of [Oil], 
and those calculated using the earlier results of McLaughlin and Bell (1998) are not 
due to relativistic effects, which are negligibly small. A fully relativistic study of 
forbidden transitions of O II appeared recently [124] . reporting a high-density ratio 
in good agreement with values deduced from the astronomical observation of planetary 
nebulae |125i 1115] . Similarly, Storey and Zeippen |126j demonstrated the necessity of 
including the relativistic corrections to the magnetic dipole operator for an accurate 

I some important corrections to the original Ml transition current |108j appeared in [109j . 
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estimation of the Ml rates for — and — ^Pi transitions within the and 
ground configurations of carbon-hke and oxygen-hke ions, respectively. The average 
value for the ratio of observed fluxes of the [O III] AA 4959, 5007 lines in spectra of 
active galaxies and quasars \T27\ supports the theoretical improvement by Storey and 
Zeippen |126] . 

The importance of the relativistic corrections to the Ml operator was also found 
in astronomical X-ray spectroscopy. The iron K lines that appear in emission in many 
natural X-ray sources [128}, have a well-known plasma diagnostics potential. In the 
calculation of atomic data for K-vacancy states in Fe XXIV used in the spectral model- 
ing of iron K lines [|128j, the situation becomes critical for the Is2s2p metastable 
level which is shown to decay through both Ml and M2 transitions. In their work, also 
based on the Breit-Pauli approximation, the authors observed that the Ml A- values 
must be calculated with the relativistically corrected operator, the difference with the 
uncorrected version reaching five orders of magnitude. 

From a pragmatic point of view, the inclusion of the relativistic corrections to 
the Ml operator established by Drake [1] in the Breit-Pauli relativistic treatment of 
complex atoms, is currently limited to SUPERSTRUCTURE |129j and to its exten- 
sion, AUTOSTRUCTURE [HI [85], using the Slater determinant approach. Other 
well known atomic structure codes allowing a relativistic treatment in the Breit-Pauli 
approach do exist (CIVS [I30], R-matrix codes [131], MCHF [1321 |133], ATSP2K 
[134] ) but systematically adopt the non-relativistic version of the Ml transition oper- 
ator [m I135[ 1136] . These atomic structure codes use the Wigner-Fano-Racah algebra 
|137[ 1138] and its extensions [139] for evaluating the angular integration of the Hamil- 
tonian and transition operators matrix elements |140[ 1141] . The starting point of this 
algebra is to find for each operator its expression in terms of irreducible spherical ten- 
sorial operators [142] . In the present work, we derive it for the relativistic corrections 
to the Ml operator introduced by Drake [1]. 
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3. The magnetic dipole transition operator 



The magnetic dipole (Ml) decay rate (in s ^) is given by |136] 



4 

3(2J„+1) 



(1) 



where u and / refer to the upper and lower levels, respectively, u is the angular frequency 
corresponding to the transition energy AEui = hw and (J/||(5*^^^|| Ju) is the reduced 
matrix element of the magnetic dipole moment tensorial operator. The latter is built 
from the vectorial magnetic dipole moment components of 



C / 2 1 2 \ 



+ '5^ 



Ze' 
mc'^rf 



Vi A iri A Si) + 



2mc^rij 



Si + 



1 



2^2,2 A (p. A 
Tjj A [Vij A {s, + Sj)] 



1 (n A rj)rij ■ {p, + pj] 



ri Apj +r,- ^pi 



(2) 



expressed in the Gaussian (mixed) system of units, as derived originally by Drake [T] 
(/^s = 2^^ )• 

In equation ([2]), one can recognize through the first two terms (L + 2S = J + S) 
the usual form of the magnetic dipole transition operator |136j . The other terms have 
been derived by Drake [1] as the 0{a'^Z'^) corrections to the usual definition of the Ml 
operator in the Pauli approximation. Some confusion propagated in the literature due 
to embarassing misprints §| but the original version [1] should be the definitive one [144j . 



This operator, written in atomic units [fiB = «/2), splits into its one- and two-body 
components 

AT AT 

Q = 5ZQ. + 5ZQ., , (3) 



i=l 



i<j 



with 

Qi = fJ'Bi {li + "isi 



aW 2d /? 



2 i (9r? ~'~ r,- (9r,- rf 80^' 



Pi A {p^ A Si) + 2r, A (n A s,) — - ^ + -s 
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Qij = flBO^ 



Vjj A [Vij ^{si + Sj)] 1 {vj A rj)rij ■ {pi + Pj) _ 1 n A Pj + Vj A Pi 



ij ij 

where e = AEui is the transition energy in hartrees. 



(4) 



(5) 



§ The two contributions in /\pj and Vj /\pi appearing with the same (— ) sign in the last term of ([2]) 
do appear with opposite signs (—)/(+) in Eissner and Zeippen and Eissner [143] . Moreover the 
factor (1/2) for the last two terms of ([2]), is missing in ref. |143] . 
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4. Derivation of the irreducible tensorial form 



4.1. The one-body operator 

The vectors Ai and Bi defined as 

Ai=PiA {pi A Si) = pi{pi ■ Si) - Si{pi ■ Pi) , 

and 

Bi = riA {vi A Si) = Ti^Vi ■ Si) - Si{ri ■ r^) , 
appearing in the one-body operator (jlj), can be rewritten in irreducible tensor form 



2V3 

2V3 



(1) / (1) (1)' 
s\ X [pI ' xpy 
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3 


yo)" 






3 



(1) f (1) 
s\ X lp\ X p\ '] 



(1) 



(1) f (1) (1) 
s; X r; ' x r 



(6) 
(7) 

(8) 
(9) 



by applying the angular momentum theory and operator techniques |145l 11421 1146j and 
using the recoupling formula for commuting irreducible tensors |147j . From the tensorial 
form of the linear momentum p = — zV (in a.u.) ^142j 



P 



(1) 



dr 



(10) 



((1) ^^"^ (1) 

Pi ^ Pi ) is written in terms of the angular part Vj^ of the 

difi^erential operator V and the renormalized spherical harmonic operator C*-^-*: 



(1) (1) 
Pi ^Pi 



(2) 



d_l 

dri ri 



(2) 



(12) 



Defining the commutator of an irreducible tensor product as in 

the two tensors ( . x j and ( x . j appearing in the second and third 
terms of ( ITTi) are related through their commutatoi!]]] that reduces to 

V6. 



iK-(v[;\c«) = (v«xc(^: 



(2) 



(2) 



-C(2) . (13) 



(-1)1+1-2 (^^(1) ^ ^( 

The fourth term of ( ITTi) is simplified by using the key reduction formula (fT^ for the 
tensor product of spherical harmonics [142] : 

k\ k k2 



-l)«v/2^( "1 " "2 IC^'^). 
^ '000 



(14) 



II We can see the general expression of the fnf " (V[2\ C^^^)) m the book of Varshalovich and al [HT] 
(equation 113 page 497). 
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Using equations (IT51) and (HM . expression (fTTj) becomes 



(1) 
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(0) 



(2) 



(2) 



(15) 



The tensor product [p^P x j?P j appearing in ([HD is simply obtained thanks to the 
relation with the scalar product Pi ■ Pi : 
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.Pi ■ p, = —A, = +— ( — + 



(16) 
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By introducing equations f|T6l) and f|T5l) in ([8]), one gets the final tensorial expression 
for Ai. 
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The tensorial form ([9]) of Bj reduces to 
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by using r 



(1) 



3 
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Inserting expressions (fTTl) and (fT8|) into equation (jl]), one obtains, after some regrouping, 
the tensorial form of the one-body Ml transition operator: 
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4-2. The two-body operator 

For the two-body magnetic dipole transition operator ([5]), we need to investigate the 
irreducible tensorial expressions of the three following contributions: 

_ Tij A Tij A (Si + Sj) 



D 



1-3 



E 



{Vi A rj)rij ■ {pi + pj) 



Vi A Pj + Vj A Pi 



(20) 



4-2.1. Tensorial form of Ci- 



Similarly to equations (IHl) and Q, the tensorial expression of the double vector 
product appearing in Cij becomes 



C 



13 



+ 



2^3 1 

3 r3. 



(0)' 



(1) 



15 1 



3 r3. 



,(2)' 



'J ^ ' ■ V ■ y " ' ij J 

Let be Xjj and Yjj, the first and the second terms of (1^ ie., Cij = Xij + Yij. 



(21) 



Calculation of Xi 



Combining the vectorial form of the tensor product of rank zero 



) v/3 
with the well-known expression |142j 



(1) 



^ 2 



(0) 



with 



e(rj - ri) 



(22) 



(23) 



k '> 

the sum over the electron pairs of the Xy contributions becomes 



(24) 



1 si Tj < 

si > . 



(25) 
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Calculation of Yij 



Starting from r*^^^ = r G*-^-* and rij = ri — r^, one first builds the irreducible tensor 
of rank two 



(1) (1) 



(2) 



y|(r?cf' + r;Cf)-2,,r,(c«xCf 



(2) 



using f |T4l) . Taking the tensorial form of r- -^ |142] 



4 >■> ■< ^ 



(0) 



(26) 



(27) 



the Yij contribution becomes 



3 rt. — T 



> ' < 



X 



+ s^P) X a 



(2) 



(1) 



) X 



(fc) 



(0) 



(1) 



3-;e7?E(-i)'i5T(2* + i)''^ 

^ ^ k 



.(2) 



(1) 



C['^ X ) 



(0) 



(1) 



3/2 



> < k 



s 



(2)' 



(1) 



X 



d''^ X c 



(k) 



(0) 



(1) 



(2^ 



The calculation of the three terms is rather tedious, requiring numerous recouplings 
for regrouping the spin and angular factors. Expressions of these terms are given in 



Appendix A The thereby defined and calculated -Fj,- and Hij contributions are relevant 



to the first and third term of (1281) . respectively, while the second term is obtained from 
the first one by merely interchanging the [ij) indices. After some further effort in 
regrouping similar contributions, the sum over all electron pairs of fl28l) reduces to 



i<j 



x<2 



i^j k 

2k 



{2k + l){k + l) r 



2k + 3 



s^> + s)'>)x Cf xC, 



(1) , o{l)^ .X (r^ik) 

2k -1 " " 
+ ( 1 - ^) ^/3{2k + 5){k + 2 



(k) 



(2)' 



(1) 



(fc+2)y2)' 



(1) 



(29) 

The final two-body contribution Ylii<j ^ijy defined in (l20l) . are calculated according 
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to equations (1241) and 



i<j 



X 



(s<" + .|")x(cf' xCf) 



(1) 



2^2 



/ fc(2fc+l)(fc + l) r'; 
(2A; + 3)(2A;- l)r'=+i 



(1) 



x/3 



X 11-4 



(2A; + 5)(2fc + l)(A; + 2)(A; + l) r| 



2A; + 3 



£:(ri - rf) 



(1) 



(30) 



4-2.2. Tensorial form of Di 



There are different ways for getting the tensorial expression of the two-body 
contribution Dij appearing in ( l20l) . One approach is to evaluate the two subparts 
i^ij ■ (Pi + Pj) and rlf'iyVi A Vj) of the operator separately. Using equations (ITOl) . (lT4l) and 
the expression of ■, C'^^^) (see the footnoted]) , one easily obtains 



rij ■ [Pi + Pj) 



d 



I r 



'drj 



dri 



drj 



d_ 

' dri 



cP X C 



(1) 



(0) 



cr X V 



(1) 



(0) 



-zv^^(q-xv. 



.(1) 



(1) 



(0) 



From the rank one tensor associated to the vector product 



Ti A Vj = -iV2 TiTj (^&P X C 

and from fl27|) and fIB.ip . we get 



(i: 



(1) 



Tj A Tj 

^3 



■^Y.^-lf^,^/{^k + l){k + l)k (Cf X Cf ) 



(1) 



fc ■ > 

By recombining equations (!33|) and (!3T|) for building the Z>jj operator, one gets 



D 



^3 



(31) 



(32) 



(33) 



d_ 

' dri 



d_ 

'dri 



^{2k + l){k + l)k (Cf ) X Cf") 



(1) 
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X 



) X C 



(k) 



(1) 



X ( Cf ) X C]^)) 



(0) 



(1) 



-^^{2k + l){k + l)k 



X V 



(1) 



(0)' 



(1) 



' i Ty 



X 



{Cf X 0]-=')'" X (Cf X V« 



(0)' 



(34) 



Let us denote the above four contributions to Dij as Kij, Lij, Mij and Nij, respectively. 
The tensor products of the second and third terms are recoupled in equations flB.2l) and 
(IB.Sp . After some tedious work, the two first contributions to Dij reduce to 



V3 
3 



k ^> 

r^k + 2 



r> 2A; + 3 

r< A; + 2 
'^2A; + 3 ^ 



+ r,, 



r> A; — 1 
2A; - 1 



d 



dr^ 



(k) 



(1) 



r> A; — 1 
'7^2A; - 1 



(1) 



(35) 



The sum of the last two terms of fl35l) being symmetric with respect to (zj)-exchange, 
the summation over all pairs becomes 



2 A; + 2 



2 A; - 1 

"■<2^ + ''>2fc^ 



)'" 



(36) 



Remembering that ^Cj^'"-' x Cj^^^ 



(1) 



C]') X C 



(fc) 



(1) 



one realizes that the third 



and fourth terms of (IMIl are also symmetric with respect to (ij) -exchange, ie. Mij 
Nji, leading to 

i<j i^j 

The final expression for J2i<i -^ij 



V3 



i<j 



j:T.(-^)'^W{2k + l){k+l)k fr 

i^j k '> ^ 



k + 2 



+ r 



- r. 



2A; + 3 >2A;-1 



X ( cf ) X C]') 



(1) 



Irreducible tensor-form of the relativistic corrections to the Ml transition operator 16 



3 tl^V V 2k + 3 r,ri+^ 

i^j k 



3 • > 

(1) 



V3>^^ , l {2k + l){k + 2) r, r^+^ 
+ 3 V 2fc + 3 r,r^;+2 



3 ' > 

(1) 



x(cf)x(cr^xv[;))^^"") 



3 Z^Z^^ ^ V 2A; + 1 r,rt+i 



J ' > 

(1) 



x(d'+^)xfc?)xv«'^'^ 



/^'(^ + 2)(2A; + 3) r 



3 V 2A; + 1 r,r^+i 



x(Cr^^x (Cf xV[,;^) ) . (37) 

4-2.3. Tensorial form of Eij 



The two terms appearing in Eij (]20l) being symmetric with respect to exchange 
{ij), one will restrict to one of the two. Using fl23|) and the tensorial form of the vector 
product, one easily obtains 



-V2YS-,f^l^^ X ((Cf X Cf>)'°' X (Cf> X Vg;)'°') 
-V25:(-l)'y2FTlA,.|_ X ((Cf X Cf )'°' X (Cf X Cj")'") 



(1) 



(38) 

The two tensor products are recoupled using equations ( 1B.3P and ( IB.ip for getting the 
following form 

Ap. _ /(E|ME5zLii X fcf-'> X (cf. X vS^"-'^ 



r,,- ^ — • ■ V 2fc + l r^+V, 



,|^,_,./5±|ME5^^ X (Cf- X (Cf X .^'f 



k > J \ 
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k > J \ 



k V ^^"^ ' '> 



k 



When substituting k ^ k — 1 and k ^ k + 1 in the last two terms, respectively, the 

/ tp,\ (k) \ ^"'"^ 

same angular contribution iC^ x Cj 1 is obtained. After substituting k ^ k + 1 
in the first two terms and summing over all electron pairs, one gets 



i<j iT^j k ^ ^ / J J 

i^j k J > \ 



i^j k 



i^j k 
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4-2. 4- Tens orial form of the two-body Ml operator 



The two-body magnetic dipole transition operator appearing in has the following 



form 



Qij = tiBOt^ 

i<j i<j 



(41) 



Replacing J2i<jCij^ J2i<j ^ij and J2i<j ^ij by their expressions ([30]), ([37]), and 
respectively, we finally get the irreducible tensorial form of the relativistic corrections 
to the magnetic transition operator 

E = -h^^' E sir, - r,) 



i<j 



X 



r- 

ij^j k * 

) + sf ) X [cr X C) J 



(0) 



(1) 



2k{k + l) 



{2k + 3){2k-l 

^2 



's'»+sf>)x[crxC'i 



(Cf xC]")''') 



(1) 



3{2k + 5){k + 2){k + l) 



4(2A; + 3) 

ij^j k 3 > 



^ ((..'" + s<")xfcl"'xCf') 



(2) 



VHk + T) 



+2 



r< k{k + 3) 
r^V 2A; + 3 

r< TA; + 2 



2A; + 3'< ' 2A;- l'> 'V 

(1) 



(1) 



d') X ( c 



,(fc+i; 



1 



(fc) 



r> V 2fc + 3 



(1) 



-2, 



(2fc + 1)2 ""r^- 



(1) 



(k-mk+2)i2k+3) f^ik+i)^ (d^)x v[;)^^'^'^ 



(2A; + 1) 



(1) 



5. Conclusion 



Systematic comparisons between different theoretical approaches and atomic structure 
codes are often used for assessing the rehabihty of the produced atomic data 
[TiSl [80] [339] . In this line, the Muhiconfiguration Dirac-Hartree-Fock (MCDHF) 
and Muhiconfiguration Hartree-Fock-Breit-Pauli (MCHF+BP) methods have been 
compared for transition probabilities in Fe IV of astrophysical interest (79] [98]. The 
authors of this comparison [98] concluded that, although progress has been made 
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since the pioneer work of Garstang |15U] , agreement between MCHF+BP and MCDHF 
values for more transitions would be desirable. No doubt that the missing relativistic 
corrections considered in the present work should be systematically calculated within 
the first method for a definitive comparison. The tensorial form of the Ml transition 
operator derived in the present paper is the starting theoretical point for implementing 
the calculation of the relativistic corrections to the Ml transition probabilities in the 
atomic structure codes based on the Breit-Pauli approximation, such as RMATRX [131j . 
CIVS [no] or ATSP2K [131] using Fano-Racah algebra [HZl [EH [152| or modern 
techniques combining second quantization and quasispin methods in coupled tensorial 
form [T58l[T5i] . 
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Appendix A. Intermediate calculations of the two-body contribution 

The rank one tensors defined by 



3 r? 



X 



It 



,(1) 



X a 



(2) 



(1) 



(0) 



(1) 



(A.l) 



<-+^v/T5 ^^^^ 



^2 /yt 2 

> < k 



=3 ' 



is; ' + sr) X (Cf^ X C, 



ft 



X 



d'^ X c 



(k) 



(0) 



(1) 



(A.2) 



and appearing as the first and third terms of ( l28l) . respectively, are transformed by 
first decoupling their spin and space parts and by using for the latter, the reduction 



formulae of tensor products involving four irreducible tensors given in [Appendix B 
These contributions can then be rewritten as: 



7(1) 



^ 3 / A;(A;- l)(2fc-3) 



10 V (2A; + 1)(2)(2A;- 1) 
A;(A; + 1) 



(Cf X Cf -^)) 



(2) 



(1) 



(1) 



[cf X +^)) 



3 / (fc + l)(fc + 2)(2fc + 5) / 
+ ''T0V (2A; + 3)(2fc+l)2 ' 



(2) 



(1) 



(A.3) 



(1) 



(2A; + 1)4 



/ (2A: + 3)(A: + l)A:(2fc-l) / (1) (1) / (,+1) 



(2)X (1) 



V5 / (2fc + 3)(fc + l)fc(2fc-l) ^ ^(1)^ ^ ^ (2)N 

5 y (2/;; -|- 1)^ y 3 \ 3 



1 /(2A; + 3)(2A; + 5)(A;+ l)(A; + 2) 



30 



(2A; + 1)^ 



(1) 



(A.4) 
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Appendix B. Reduction of tensor products involving four irreducible 
tensors 



((Cf xCfl)'°'x(cf' xCf)"') 



(1) 



V6 I {2k- i){k - i)k f^^k-i) ^ 



6 V (2fc + l)3 V ^ 



VQ I {2k + 3) (A; + 2) {k + 1) (^^k+i) ^ 



6 V {2k + lf 



\/3 /(2/c-l)(fc-l)(/. + l) 



3 V {2k + lf 

'cf+^^xCf+^M^' (B.2) 



x/3 /(2A; + 3)(A; + 2)A; ^^(,+1) 



((Cf)xcf)^^(cf)xvg))^°^) 



(2A; + 1)2 

(1) 



/ (2A:-l)(fc-l)(fc + l) / / (,) -(DXe^-i)^ « 



(1) 



V3 / (2fc + 3)(A: + 2)/c / / (,) (dxC^+dV^) 

In the first two formulaes, the tensorial operators Cf'\ Cf''' et Cf \ Cf^ act in 



different spaces and all commute with each other. To get expressions ( IB.ip and (1B.2P 
we have used the following transformation [147] : 



Cf^xcf)^'^^ 



, , (1) { k k k^\ 

(cf^xCf) =5^(-irV5(2^7 + l)(2/. + l)<^ 1 1 A:2 [ 

^ [ (7 /i 1 J 



xffcf)xCf))^^^(cf)xCf)^^^V' (B.4) 
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where the cases {ki = 0/k2 = 1) and {ki = l/k2 = 0) correspond to (IB.ip and (lB.2p . 
respectively. The sum over g and h is hmited to g = h = \k — 1\ and g = h = k + 1. 
By applying (HM to both tensorial products appearing in ( 1B.4I) . we get the right-hand 
side of equations (1B.1|) and (IB. 21) . To obtain (IB.Sp . we also use the transformation (IB. 40 



replacing Cj^^ by the operator Vq'' , keeping in the sum the {h = k — l/g = k — 1, k) 
and {h = k + 1/g = k, k + 1) contributions. 
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